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Biorthogonal Laurent polynomials, Toplitz determinants, 
minimal Toda orbits and isomonodromic tau functions 
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Abstract 

We consider the class of biorthogonal polynomials that are used to solve the inverse spectral problem 
associated to elementary co-adjoint orbits of the Borel group of upper triangular matrices; these orbits 
are the phase space of generalized integrable lattices of Toda type. Such polynomials naturally interpo- 
late between the theory of orthogonal polynomials on the line and orthogonal polynomials on the unit 
circle and tie together the theory of Toda, relativistic Toda, Ablowitz-Ladik and Volterra lattices. We 
establish corresponding Christoffcl-Darboux formulae. For all these classes of polynomials a 2 x 2 system 
of Diffcrential-Difference-Deformation equations is analyzed in the most general setting of pseudo mea- 
sures with arbitrary rational logarithmic derivative. They provide particular classes of isomonodromic 
deformations of rational connections on the Riemann sphere. The corresponding isomonodromic tau 
function is explicitly related to the shifted Toplitz determinants of the moments of the pseudo-measure. 
In particular the results imply that any (shifted) Toplitz (Hankel) determinant of a symbol (measure) 
with arbitrary rational logarithmic derivative is an isomonodromic tau function. 
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1 Introduction 

The connection between orthogonal polynomials on the line and Toda lattices is rather well known [3] , as 
well as the relations to the KP hierarchy . Dynamical variables of the Toda lattice are arranged into a 
tri-diagonal Lax matrix, that can be viewed as a recurrence matrix for a system of orthogonal polynomials. 
In the (semi)finite case, the evolution of the corresponding measure provides a linearization of the Toda 
flows. More generally, one can set-up (in)finite-dimensional Hamiltonian systems on R 2 ™ (n < 00) with 
Hamiltonians 

1 n \I\ 

i—l i^I j=l 

I:={h <i 2 < ... <i k } ■ (1-2) 

As it is noted in JU] such family of Hamiltonians (labeled by the multi-index I) contains integrable lattice 
hierarchies of Toda, relativistic Toda, Volterra and Ablowitz-Ladik type. These integrable Hamiltonian 
systems have a Lax representation with Lax operator given as a n x n lower Hessenberg matrix which we 
denote by Q (in |5J^] it was denoted by X), belonging to a certain "elementary" (2n — 2)-dimensional 
co-adjoint orbit of the solvable group of upper triangular matrices. These systems are linearized by the 
Moser map 

Q - W(z; Q) := (si - Q)^ = £ Mgl (1 _ 3 ) 

3=0 

Aj(Q) = Q j n (i-4) 

In the case of infinite lattices these expressions take on a formal meaning in terms of power series but 
the analysis is unchanged. 

The moments fij of Q define a normalized moment functional C and the reconstruction of Q from its 
moments (the "inverse moment problem") can be accomplished by constructing a suitable sequence of 
biorthogonal (Laurent) polynomials {ri,pi}i e N 

£(npj) = , 

where p^s are polynomials in x of degree i while r-j's are, in general, polynomials in x and a; -1 . The 
(infinite) Lax operator Q corresponding to the chosen orbit is then reconstructed by ^U] 

Qij = C(rixp 3 ) . 



2 



Explicit formulas for these biorthogonal polynomials in terms of shifted Toplitz determinants can be 
found in (91 110) and will be recalled here in due time. Vice versa, one could assign an arbitrary (generic) 
moment functional C : C[z,z _1 ] — > C, a multi-index / and then reconstruct the Lax operator Qi (i.e. 
view the Lax operator as a function of C rather than the other way around) 

C^Qi(C). 

From this point of view, the linearization of the (infinite) Hamiltonian hierarchy is accomplished simply 
by 

£ t («)= J C(eSiV*M > .) ) (1.5) 

where the series may have to be understood formally. This procedure displays the common nature of all 
the above-mentioned integrable lattices, inasmuch as the linearizing space is always the same (the space 
of moment functionals) and what changes from one lattice to another is only the orbit, namely the map 
Qi- 

Finite dimensional systems(of dimension 2n — 2) on an elementary orbit Qi correspond to those 
moment-functionals for which certain shifted Toplitz determinants of size < n do not vanish whereas all 
larger ones do. In such cases, the tau function of the hierarchy defined by the (closed) differential 

n , 

dlnT = ^-Tr n (Q 7 )dtj (1-6) 

.7=1 J 

and coincides with the largest non- vanishing (shifted) Toplitz determinant. 

One of the main purposes of this paper is to connect this determinant to a different notion of "tau" 
function, namely the one introduced by Jimbo, Miwa and Ueno in [131 114) . It was shown in (S) 0] 
that the Hankel determinants of an arbitrary (generic) "semiclassical" moment functional on the space of 
polynomials can be identified with the isomonodromic tau function introduced by our Japanese colleagues. 
Similarly, it was shown in |16) that Toplitz determinants of a particular class of symbols on the unit circle 
are also identifiable with the same kind of isomonodromic tau functions. 

These two apparently distinct situations are in fact the two ends of a "continuous" spectrum of 
situations: in fact the case of Hankel determinants is dealt with in the setting of (generalized) ordinary 
orthogonal polynomials, whereas that of Toplitz determinants uses orthogonal polynomials on the unit 
circle; in this latter situation one considers polynomials Pi(z) orthogonal in the usual L 2 (S 1 , d/i) sense 



Pj(z)p k (z)d(J,(z) = 5jk ■ (1-7) 

s 1 

Here one defines rj(z) = pj(z _1 ) and the orthogonality is recast into 

C(r jPk ) = S jk , (1-8) 

where - in this special case - 

C : C[z, z- 1 } -» C ; C(z 3 ) = [ z J d/i(z) . (1-9) 



We see that we can regard the case of orthogonal polynomials on the circle as a special case of biorthogonal 
Laurent polynomials with respect to a moment functional satisfying the reality condition fi k = M-fc- 

According to the previous description of integrable lattices, the two situations correspond to two 
different elementary orbits and hence we should be able to treat them on a common ground, together 
with all the other lattices associated with the orbits 0/. Indeed, we will show that this is the case and 
that for the class of moment functionals of the semiclassical type introduced in ^ all the shifted Toplitz 
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determinants which arise as tau functions of the corresponding integrable lattices are also isomonodromic 
tau functions for a rational 2x2 connection on C 1 which will be explicitly constructed in the paper. 

The approach to this problem follows the strategy used in 0] rather the one in JHj; m the course of 
our analysis we will obtain generalized Christoffel-Darboux identities which naturally interpolate between 
the ordinary CD identity for orthogonal polynomials on the line and the one for orthogonal polynomials 
on the unit circle. 

Moreover we will show that the Toplitz and Hankel determinants of the same size for one such moment 
functional are connected by a sequence of elementary Schlesinger transformations, at each step of which 
we obtain tau functions associated to interpolating orbits; in figurative terms, we show that the papers 
[TB| (see Example 110. Q and jjj are connected by a Schlesinger transformation (when specializing the 
semiclassical measure to the one relevant for ^B]) and that "neighboring" elementary co-adjoint orbits 
are also connected by an elementary Schlesinger transformation. 

Acknowledgements. The authors thank John Harnad for stimulating discussions. M. G. is grateful 
to Laboratoire de Physique Mathcmatique, Centre de Recherches Mathematique and the Concordia 
University for their hospitality during his visit to Montreal, where the work on this project has started. 
During the later stages of preparation of the manuscript he also enjoyed hospitality of Institut des Hautcs 
Etudes Scientifiques. 



2 Setting 

We start in the most general and abstract setting, without any reference to a (pseudo) measure. We 
consider an arbitrary moment functional 



C : C^z" 1 ] 



(2-1) 



on the space polynomials in z and z 1 and denote its moments with fij = C(z 3 ) , j € 1. We introduce 
the following shifted Toplitz determinants and polynomials 



A* 



det 



fit-i 



fit 



A* 



fl£+n-l \ 
fli+n-2 

fit 
Ai E 



(2-2) 







fii-i 



p e n (x) := det 



fit 



fli-n+l fl'i-n+l 
V 1 X 



fi>e+n \ 

flt+n-1 



M-e+i 
x n 



(2-3) 



Using some classical identities for determinants we can derive recurrence relations for the shifts n — > 
n + 1 and I — > I + 1 for the above polynomials. We first need the following 

Proposition 2.1 For any (n+1) x (n+1) matrix A the following determinant identity holds true (Jacobi 
identity) 

il..n i2..n+l a2,,ti+1 a 1..TJ A 1. .n+1 a 2. .n jr. A \ 

/i l..n /i 2..n+l ~ A l..n ^2. .n+1 ~ /1 l..n+l /i 2..n 7 l z_4 i 

where the sub/super-script ranges denote the rows/columns of the submatrix we are computing the 
determinant of. As a corollary, for any (n + 1) X (n + 2) matrix B we have 



pl-n+l dL.ii 
D 2..n+2 D l..n 



,l..n+l 



,B 9 



r>L..n+L Tjl-.n _ pX. 

JD l..n+l r '2..n,n+2 — - D l..n,n+2- u 2..n+l 

which can be obtained from (|2-4|l by adjoining an appropriate row. 



(2-5) 
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Using 1)2-4(1 on the determinant defining p l n we find 

X&nPn-l - A n Pn-1 = A n-lPn 



(2-6) 



Applying (|2-5J) to the determinant defining adjoined of the next row of moments on the top we find 



, 1-1 \i i \L , e-i , i \t-i 

Pn A n + A n+lP n -l = Pn A n 
Pn A n +xA n+1 p n _ 1 = p n A r 



X A n Pn-l ~ A n Pn-1 = Pn A n-l 

We now use these identities to express p l n := [p l n , p^-i] m terms of p t n _ 1 — [pl_u P^X] 

xn -n ^n ^n-1 ^„ 



(2-7) 
(2-8) 
(2-9) 



Pn 
Pn-1 



I _ nl I 
Pn ~ L nPn-l 



AU (Af^) 2 (A^r) 2 

At\ A» 

A 1 i A £ , 

n— 1 n— 1 

£ ^2 



Pn-1 
Pn-2 



dctCf, = 



(A„) 



(K-i) 2 
]Ci(x)- x Ci{y) -j=(l 



Circle Case 



A 



ra-1 







A, 
-1 



where 



J := 



1 

-1 



(2-10) 

(2-H) 
(2-12) 

(2-13) 
(2-14) 



We have named this the "circle case" because this sort of recursion is relevant for orthogonal polynomials 
on the unit circle. We next derive a recursion in I 



Pn-1 



A 



A 1 A e_1 A e A 

^n+l^n-l ^n+l^n 



A 



xiK-y x(Al" 1 )2 
A e -\ A e 



Pt 1 

Pn-1 



(2-15) 



t _ rrl t-1 
Pn 1 n Pn 

1 (Al) 2 



detTf 



x (Al" 1 ) 2 



jr„ f (x)- 1 r„ £ (y)-j= i 



A 



n-l 



1 



A! 




(2-16) 

Circle to Line Transform (2-17) 

(2-18) 



The name "circle-to-line" refers to the fact that this recursion relation interpolates between the previous 
"circle" case and the next one, which will be named the "line" case. Indeed, composing these two we can 
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express p l n = [pi, p K n _\] in terms of p e n _\ = [p l n _\, p l n l%] 



-1 , 1-2 - 



Pn 
Pn-1 



n — 1 



(Al) 2 , 1 

fA* -1 ^ Pn-1 
L Pn-2 



Line case 



dot £ 



(Af-\) 2 
]Ci{x)- x Ci{y)-i = {x-y) 



A 







o o 



(2-19) 

(2-20) 
(2-21) 

(2-22) 



This recursion is called "line" case because it is the relevant recursion relation for ordinary orthogonal 
polynomials on the line. 

2.1 Second— kind polynomials 



Let us define the following second-kind polynomials 



<{x) = C z 



PnO) - pi(z) 



(2-23) 



The three types of recursion (|2-10l 12-151 12-19J) involve at most a multiplication or division by x and have 
otherwise constant coefficients (in x): moreover we find 



x-'Kix) 



x — z 

r - Z ~ X P 1 A Z ) \ 1 r i -l e f „ss 

£zl ) - -C z {z p n (z)) 

x — z x 



(2-24) 
(2-25) 



The last terms in these identities vanish because of the determinant structure of p n , provided that n > 1 
and < £ < n — 1 for the first case and — 1 < £ < n — 2 for the second case. From this observation we 
find that these auxiliary sequences of polynomials satisfy the same recurrence relations in the following 
ranges 



K 
nt 



= < L n 



n-l 

-2 
n— 2 



n-l 
-1 
n—1 
1 

n-l 
2 

n-2 



l<£<n-l 
0<£<n-2 
0<£<n-l 



(2-26) 



3 Christoffel-Darboux formulae 

Consider (n,I)eNxR and choose an arbitrary path starting at the origin of the following type 

{(n fc , 4), * = 0, 1, ... , (n , 4) = (0, 0), (m, £i) = (1, 0)} (3-1) 
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and such that the possible subsequent moves are right, up or up-right. For the move (n k -i,£ k -i) > (n k £k) 
we introduce the transfer matrices following an idea of ^2] used for orthogonal polynomials on the 
circle 

{C„k if the move is right (circle move) 
if the move is up (circle-to-line move) (3-2) 
L ik k if the move is up-right (line move) 

Using these transfer matrices we define the two dual auxiliary sequences of matrices as follows 



E k (x) = T k (x)^ k -i{x) 



detT fc (z) 



"0 — "0- 



(3-3) 
(3-4) 

(3-5) 



This definition in particular implies that 



n.trdetr, 



(3-6) 



The choice of the initial conditions for the auxiliary sequences is arbitrary but it is convenient to choose 
S in such a way that the first column of E n will contain <p l n and p^-i an ^ the secon d column the 
corresponding second kind polynomials. Since the matrices constructed with the polynomials p * and 
the second kind polynomials already satisfy the same recursion relation for k > 1, it is sufficient to impose 
the same initial conditions with the following choice (recall that the first move is always a circle-move) 



"0 



n 



m — 1 



i 



1 /Hi - hqx 



fi x - fix /i 
1 



(3-7) 



Recall that for any 2x2 matrix we have A f — det(A)jA x j 1 . We now compute 

~t(x)j~ k (y) = 1 r Ei_ 1 (x)T t k (x)jT k (y)E k _ 1 (y) = 
det Jfc(x) 

= Zt_ 1 (x)jT- 1 (x)T k (y)Z k _ 1 (y) = 

= Ht^jjjBufel + StiW (jT^(x)T k (y)-i) E k ^(y) 



(3- 



Let us define l k := l k 
uniformly written 



l k -\ and h k := n k — n k -\. Then the three formulae I|2-1UI2-15I2-19|) can be 



jT- 1 (x)T fc (y)-j = (-l) 1 



x tk y nk 



n k -l 

4-1 



1 - h k 



(3-9) 



detT fc (x) = {-l) n «x 



Air 1 



Y[detT 3 (x) = (-ir k x n "- e " (A& 

3=1 



(3-10) 
(3-11) 



7 



Summing up both sides of eq. (|3-8|1 we obtain the following master Christoffel-Darboux identity 
E* N (x)jS N (y) 



-1/y 
1/x 



y x 



l l 

x y 



N-l 



k=0 



N-l 



k=0 



A tfc + l-J- 

A^fc+1 
^Hfc + l 

4+i - l 

A 4+1 

4+i - 1 



1 - fik+i 



l - ^fc+i 



(3-12) 



3.1 Principal CDI 

We look at the (1,1) entry of the above identity 

<-W.fc--teU w _^W 



(Al™ ) 2 I x nN ~^ N bnNyl " x 7lN ~ 



1 1 

a; y 



JV-1 

E 

fe=0 



(A^) 2 



A«fc+1--L 

A £fc+1 



pt(y) + ( 1 -^+i)^ f l- 1 i(y) 



P&O*) - (i - 4+0-j^tp^- 1 !^) 

^n fe+ i-l 



(3-13) 



The two terms in the product inside the sum here above can be simplified using (|2-7[) for the case 4+i = 
and (|2-9|l for the case n^+i = indeed 



j^pt(y) + (l-n k+1 ) P i«-\(y) 



^lP P n k k+1 -M ifn fe +i=0 



= < 



n fc + i-l 



, P&(l/) ifn fe +i = l 

.e k+1 Pn k+1 -i\y) 
^"k+l 



(3-14) 



A +1 

P ^(x)-(l-4+i)-^ r p^ 1 1 W 
^nfc+i— i 



A fc+1 
A^+i- 

"n it +i — l 



rP«t +1 (*) if 4+i=0 



Pi\(x) = p e n \ +1 - 1 (x) if 4+i = l 



A 



1 P"fc 



A £fc 



Using these expression in the RHS of (|3-13(1 the identity becomes 



(-i)" w [pftpiOs) 



p&(y) 



r n N -l N b "ri N 



PnV-M 



N-l 



(l-OEf-i)"' 



fc=0 



A 4+1 A 

^"fs + l^nfc 



(3-15) 



(3-16) 



We can repeat the same arguments for the second-kind polynomials appearing in the other matrix entries; 
care must be paid to the fact that (So)i2 is not TZq = 0. 

We obtain the following supplementary CDI's (provided that < ik < rik+x — 2, k = 1, . . .) 



(-i)" w ( <V-\{*) 

(Al~) 2 



pn N (y) 



r n N -£ N 5 "jv 



pf-i(v)]-z = 



N-l 



(f-i)£(-D- 



fc=0 



A^fe+l A t k 



\2 



(I" 1 ) 



'V-i (v\in k+1 ~ (t\t^+^ 



\2 



JV-1 



fc=0 



x n N -i N '^nM-lKV) J - 

A^fe+l A£fc 



(3-17) 



The additional term in the second identity stems from the mentioned discrepancy in the definition of so 
with the definition of the auxiliary polynomials: indeed the term with k = in the sum (|3-12|1 is not zero 

Thus the second identity above is rewritten as 



in the off-diagonal terms but 







\x /A 1 ^ A 4+i A 4 



(3-18) 



4 CDIs for biorthogonal Laurent polynomials 

The formula? derived in the previous sections for the Christoffel-Darboux identities are very general 
however the (Laurent) polynomials that appear in the sum are not biorthogonal with respect to the 
moment functional C unless the sequence rik is strictly increasing and the sequence £k is weakly increasing. 
This is the situation which interests us the most and hence from now on we will assume that rift = fc 5 . 
Moreover all the elementary orbits of the integrable lattices we are considering are in correspondence 
with this situation. 

From the formulae defining the polynomials p l n it follows that 

A (p£(*)p£ +1 -W m+l ~ m ) = W-irAfrA&i 1 ■ (4-1) 
This suggests that we introduce the following monic polynomials 

7r„(s) = -j-ffcix) (4-2) 
A„ n 

5 If n k were not strictly increasing then the polynomials would be biorthogonal only provided the moments satisfy some 
non generic condition of vanishing of certain determinants. 
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(— 1 Y n 

Pn ( x) = y_^i-» p ^-i (it) . ( 4-3) 

An 

It is understood that the determinants A„™ must not vanish: this is our implicit assumption of genericity 
on the moment functional. While the 7r„'s are monic in the usual sense, the p n : s are normalized on cither 
the highest or the lowest power depending on l n +i- Moreover the 7r„'s are polynomials in x whereas the 
p n 's are polynomials in x and x~ x . They satisfy the orthogonality relations 

A £ " + i 

£ z {Pm{z)-K n {z)) = 8 mn h n , h n := " + . (4-4) 

A„" 

We finally introduce the (bi)-orthonormal polynomials and the second kind polynomials 

l $h 

p n (x) := —=n n (x) = " 

V"n ,/a'«A" +1 



p n (x) := C z 



Pn{x) -Pn{z) 



x — z 



r n (x) := -±= Pn {x) = Y P \' X 



and their "starred" 



Pl{x) := x^- n+1 p n (x) p* n {x) := ^"-" +1 p„(x) (4-6) 
<(x) := x n - e ^r n (x) r*{x) := x"- £ "^f„(x) (4-7) 

In terms of these (Laurent)polynomials the CDIs read 

JV-l 

r n{x)Pn{y) = IN (pN(y)r N -l( x ) - P*N(x) r *N-l(y)) 

n=0 
N-l 

(V- x ) r n (x)p n (y) = (pN(y)r N -i( x ) -PN( x ) r N-i(y)) + 1 

n=0 
N-l 

{y-x) Y r n {x)p n (y) = IN (p N {y)r N -i(x) - P N ( x )rN-i(y)) - 1 

n=0 
JV-l 

{y-x)^2r n (x)p n (y) = 7jv (p N {y)r N -i(x) - P N (x)r* N _ 1 (y)) 



n=0 



It is convenient to rewrite in matrix form the previous identities as follows 

p(a;) := bo, ...]*, p(a;) := [po, ■••]*, r(x) := [r , ■■■]*> r(x) := [f , ■ ■ ■]' (4-9) 
P(x):=[p(x),p(x)} , R(ar) := [r(a;),f (a;)] (4-10) 
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N-1 

(njv-i).y := V S ik 8 kj (4-11) 



k=0 



R £ (x)n JV _ 1 P(y) = lN 



p* N (x) rrf-i(x) 
P* N (x) rN-i(x) 



Pn{v) Pn(v) 

r *N-i(y) r *N-i(y) 



(4-12) 



Remark 4.1 A word about the relations with previously known (hi) -orthogonal polynomials is now in 
order. If all the moves (except the first one) are line- moves namely if £ n = n — 1 then it is not hard 
to show that ir n = p n are just orthogonal polynomials with respect to the (restriction of the) moment 
functional C to positive moments. Moreover the shifted Toplitz determinants A"~ x are (up to a sign) the 
same as the Hankel determinants of the same size (by permuting appropriately the columns). 

Vice-versa, if all moves are circle-moves (i.e. l n = 0) (and we also impose certain reality conditions 
on the moments of the functional) then the 7r n are orthogonal polynomials for a certain measure on the 
unit circle and the p n 's are their so-called "dual" Laurent polynomials. The determinants appearing then 
in our sequence are precisely the "standar" ones A° . 

A second remark is that all these polynomial do satisfy three-terms recurrence relations, although of 
a different sort than the standard ones. Indeed, it is well known that orthogonal polynomials p n satisfy 
relations of the form 

Xp n = InPn+l + PnPn + Jn-lPn-1 , ( 4 "13) 

where the coefficients 7„,/3 n enter in the tridiagonal Jacobi matrix representing the multiplication by x in 
the basis of the p n 's. At the opposite "end of the spectrum" , orthogonal polynomials on the circle satisfy 
a different sort of three term recurrence relation, of the form 

x(j>„ + 5nPn-l) = InPn+l + PnPn ■ (4-14) 

It is not hard to show J!A \l(Jf that the polynomials that we are considering precisely "interpolate" these 
two sorts of recurrence relations as follows 

X(p n + (1 ~ in)S n Pn-l) = InPn+l + inPnPn , ( 4 "15) 

for certain coefficients "fmPn,S n whose explicit expression in terms of Toplitz determinants can be ob- 
tained from the formulas above but is irrelevant for this discussion. We see that "circle moves" (£ n = 0) 
correspond to a three-term recurrence relation of the type appearing for O.P. on the circle, while "line 
moves" (£ n = 1) correspond to the "usual" recurrence relation. 



5 Infinitesimal deformations of the moment functional 

We study the infinitesimal deformations for the wave vectors p(x), p(x), r(x) andr(i) under an infinites- 
imal deformation of the moment functional. Let us introduce the matrix of recurrence for these sequences 
of polynomials 

xp = Qp ; a;r* = r*Q , Q nrn := £(zp n r rn ) . (5-1) 

The matrix Q is of Hessenberg form, namely has nonzero entries on the superdiagonal and possibly on 
the diagonal and all other nonzero entries in the lower triangular part. The biorthogonality relation can 
be rewritten as 

C [pr ( ] = 1 (5-2) 
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Suppose we infinitesimally deform the moment functional 



£(.) = -C(F(z) 



(5-3) 



Here F(z) can be any function (even a generalized distribution as we will see) provided that the moments 
of the deformation are still well defined: if C is given by an analytical expression in terms of some integral 
representation (as we will assume later on) then this means some condition of analyticity on F: if the 
functional is only defined by its moments, then F should be interpreted as formal series. In any situation 
the typical case of F being a polynomial (corresponding to the usual formal Toda-type flows) will be well 
defined. 

A little more generally we could even assume that F is a distribution, particularly delta functions or 
derivatives of it. For instance we can consider deformation of the type 



5C( P (x)) = C(p(x)) = - 



d 

da; 



p(x) 



(5-4) 



for some constant a: this means that we (formally) have set F to be the fc-th derivative of the Dirac delta 
distribution for the given moment functional supported at x = a. 
Corresponding to any of these deformations the BOPs deform as 



tfp = U (F) p , Sr = U (F) r 



(5-5) 



where a priori U and U are lower triangular matrices since the range of powers of x entering in the 
expressions p n , r n will not change. In order to find expressions for these matrices we note first that their 
diagonals are the same 

1 



(U (F) )nn = (U< 



= --Sln(h n ) 



Indeed we have 



Sp n = S—= + 



= — — 5la(h n )p n + previous 
1 



Sr n = ——51n(h n )r n + previous . 

Differentiating the orthogonality relation we obtain 

F(Q) for the case of an ordinary function F 



(5-6) 

(5-7) 
(5-8) 



) p(a;)r*(aO 



(5-9) 



for a deformation supported at one point 



and hence according to the two types the matrices describing the infinitesimal deformations are given by 

Hi?) = F(g)_ , U (F) = F(QY_ (5-10) 

U«) = d k a (p(a)r*(o))_ , U^> = d k a (r(a)p*(a))_ (5-11) 

where A_o means the lower triangular part plus half of the diagonal. Note that from l|5-()l) and the 
definition of h n it follows 



k=0 
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<5/lnA^ =-Tr„F(Q) 



5s* In A£ = -d*J2Pi( a ) r i( a ) 



(5-12) 



where we have used the notation for the truncated trace Tr n A := X)j=o ^ 



33 ■ 



5.1 Deformations for the second-kind (Laurent) polynomials 

Using Leibnitz's rule we obtain the following deformation equations for the second-kind wave vectors p, 
r. For a deformation by a function F{x) we have 

s F , - K) _, w)g+£s (M) p _(M)) ei 

while for F = 5^ (z — a) we have 



x — z 

F(x)-F(z)\ r _ ( F(x)-F(Q*y ]c 



x — z 



(5-13) 



6 F p = U^p- 



d k p(x) — p(a) 
da k x — a 



s F r = m*- jl£Wz£W 
6 Folded version of the deformation equations 



(5-14) 



Let us define 



Xn 



Pn Pn 



(6-1) 

' ri-1 'n-1. 

We want to express the previous infinite-dimensional deformation equations in terms of Xn alone; this 
process is conceptually identical to the one followed in [3] and which is named " folding" . To this end we 
formulate the following 

Theorem 6.1 The infinite deformations $5-10\) for the wave vectors p, r and for the second-kind wave 
vectors p, r 15-131 \5-14\j are equivalent to the following deformation equations for x n , n > 1. 

5 { F)Xn = Ui F \x) Xn+Xn U (F) ' R {x) 

5 {st)X n{a) = Up(x) X n(x)+Xn(x)U^' R (x) 
where we have used the following definitions: 



(6-2) 



\F{Q) nn 

F(x) - iF(Q) n _i in _! 



U (F),R 



W F 
-F(x) 



-In 

F(x) - F(z) 
x — z 



-(^qF) (n-1)*, n-1 (VQf)n,n-l 



Fix) - F(Q) , , 
V Q F := (6-3) 



Up{x) 



gk 1 
daJ2 



PnTn 

o 





-Pn-Vn-1 



da k x 



-Pnfn-l PnP n 

-fn-iTjt-i r n -ip n 
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u^ R (x) = d k a 







(6-4) 



Here, for a function f(z) we have set 

f(Q)i^ := £(nf(z)r*) , f{Q)i*j := C{p* f(z) Pj ) . 



(6-5) 



Proof. We compute the deformations of both rows of Xn- We start with deformation involving a function 
F(x): the first row deforms according to the equation 



S F \p n (x),p n (x)} = 5 F e l n • [p,p] = 4 • U (F) • [p, p] + e l n • [p.p] 
where we have set 



W F 

-F{x) 



x - Q 



W F (x) :=C z 



F(x) - F{z) 



x — z 



(6-7) 



We now note that 



4 -U^ • [p,p] = ^F(Q) nn [ Pn ,p n ] +e l n C z (F(z)p(z)r t (z)n n _ 1 [p(a;),p( a; )]) 
= \F{Q) nn \p n ,Pn] + e* n C z ((F(z) -F(x))p(z)r*(z)n n _ 1 [p( 2 ;),p(x)]) 



= ^F(Q) nn [p n ,p n ]+e t n £ 

~^-^{Q)nn[Pn-i Pn] ~~ 'Jn^z 

This implies that 



F(z)-F(x)_ 



x — z 
F(z) - F{x) 



P{z)(ln\Pn{z)*,r„-l{z)]jXn{x)- [0,1] J J = 

Pnll'n-r,. l ])iXn(x)+e t n . F{X ll^ Q) -[0, ei ] (6-8) 



Sf \Pn(x), Pn (x)] = 2 F (Q)nn[Pn,Pn] ~ In^z 



F(z) - F(x) 



z — X 



Pn\Pn,r n -l] I JXn(x) + 



+ e„ • p.p 



W F 
-F(x) 



(6-9) 



In a similar way we can compute the following deformations 
6 F [r„_i(x),f„_i(x)] = Spe^ ■ [r,r] = 
= et 1 -U^-[r,7] + •[!-,¥] 



W F 

-F(x) 



-n-l 



x-Q* 



[0,ei]. (6-10) 



The computation now involves 



eti-U< F Mr,rl 



F(Q)^i ir ^i[r n _i,r„_i] + • F(Q*)n n _i[r,r] = 

,n— 1 [^*n— 1 j ^n— 1 ] + e «-i£z (^( 2; ) r ( 2; )P*( 2; )nn-i[r(a;),r(a;)]) 



= (f(x) - ^F(Q) n - hn -^ [r n _i,r„_i] + 

+ eti£, ((F(z) - F(x))r(z)p*(z)n„_ 1 [r(x),?(x)]) 

= (f(x) - ^(Q)„-i,„_i) [r n _i,F„_i] + 
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F(z) - F(x) 



^n—l^-'z 



F ( x ) - ^F{Q)n-i,n-i ) [r„_i,r„_i] + 



SF(z)-F(x)_ 



r(z) - 7»M*).<-i(*)]jx£(aO + [0, 1] 



F(x) - F{Q t ) 



z — x 



r„_i[p„, r*^] ) $X*n( x ) + e «-i 



ei 



where we have used the following definition 



An ^ An 



Thus we have obtained the following deformation equation 

5 F [r„_i(x),f„_i(a;)] = (f(x) - ]^F(Q) ri -i. n ^ [rv_i,r„_i] + e*^ • [r,r] 



W F 
-F{x) 



F(z) - F{x) 



z — x 



JV-l[Pn,<_j] )jXn ( X ) 



By "starifying" both sides we obtain 



Wf 

-F(x) 



( F(z) ~ F{x) 



\ z — x 

Putting together l|t)-9|) and (|(j-14|) we obtain finally 



r n -i[Pn-,r* l _ 1 ] JXn(x) 



$FXn 



Ui F Xx) X n+XnU {F) ' R {x) 
\F{Q) nn 







F{x) - ^(Q)^!,^! 



\F{Q) nn 

F{x) - ^(Q)^!,^! 



In 



F(x) -F(z) 



PnPn PnTn-1 
Pnfn—1 f n —l r n—l 



— (VqF) u ^i (VqF)„,„* 
-{^QF)(n-l)*,n-l (Vgf)n,n-1 



W F 
-F(x) 



We now consider a deformation supported at one point z = a. 



S F \p n (x),p n (x)} = ei ■ U^) • [p,p] - e*„ ( - 



p(x) - p(z) 



x — z 



[0,1]. 



This time we have 



e„ • U<<> • [ P ,p] = i^(p„(a)r n (a))[p n ,p„] + ■ p(a)r t ( a )n„_ 1 [p, p] = 



1 



9a(j>n(a)r n (a))\p n ,p n } + d a e n ■ p(a 



7»bn( a )^n-i(a)] .. r x [0, 1] 

J An I, X ) 



We thus have 



Sf [p n (x),pn(x)} 



rO a (p„(a)r n (a))[p„,j) n ] + <9 a JXn{x) - d a 

Zi X CL X Qi 
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Similarly for the Laurent polynomials 



S F [^(aO.TVxO*)] = • U<0 ■ [r,r] - ( ^ 



r(x) — r(z) 



[0,1] 



where now 



-<9£(Pn-i(a)r^i(a))[r^i,r n _i] + 0£<_i • r(o) 



et 1 -U(^.[r,?] = 

7n[Pn (a),<-i (a)] . *t, \ [0,1] 
J X n W 



so that finally 



(6-22) 



(6-23) 



Sf [r n -i (x) , r„_i (x)] = - - d% {p n -i (a)r„_i (a)) [r n _i , r„_i] + 



+9, 



■ e 7n bv-i (a)p„ (a) , r„_i (o )r^_ x (a)] . ^ t _ , [0, r„_i (x)] 



x — a x — a 

Starifying this last identity and collecting it together with l|6-21|) we finally have 



d k fir Pn(a)r n (a) 

-p„_i(a)r„_i(a) 



<9a fe I 2 
d k 1 
d^2 





Pn^n 

-_p„_ir„_i 



Pn(a)Pn(a) pn{a)rn-i(a) 
r n _i(a)p„(a) r„_i(a)r^_ 1 (a) 



3* 7n 

9a fc x — a 



-i- 

a— a.' 

o o 



(6-24) 
(6-25) 

M 

(6-26) 
(6-27) 



This concludes the proof. Q.E.D. 



7 Moment functionals of integral type and ODE 

We now assume that the moment functional that we are considering admits an actual integral represen- 
tation 

£( z fc ): = V^/ e- v{z) z k Az. (7-1) 

As far as the previous discussion on deformations is concerned, the integral representation of the moment 
functional is largely irrelevant, the only issue being the convergence of the deformation function: therefore 
the "potential" V(z) as well as the sets of integration Tj could be completely arbitrary. However, in view 
of our intentions, we will assume that Tj are contours in the complex plane and that V(z) is a locally 
defined smooth function on these contours with the only restriction coming from the fact that negative 
moments should be defined as well as the positive ones. 

In fact -although many considerations would remain identical in more general situations- we will 
assume that V is a locally analytic function in the complex z-plane excepted at some punctures, identically 
to the case of semiclassical moment functionals studied in [51 E] with the only extra restriction that 
all negative moments should be defined and finite. 



16 



Semiclassical Moment Functionals. For the reader's convenience we briefly recall how these semi- 
classical moment functionals are constructed In this case the potential is such that the 
derivative is an arbitrary rational function 

V'(z) = Rational function (7-2) 

and thus V(z) is a rational function plus logarithmic singularities at those poles of V where the residue 
does not vanish. For simplicity we assume that V has either a pole or a nonzero limit at z = oo. Once 
we have chosen the potential V we also choose an arbitrary collection of contours (avoiding z = 0) {Tj} 
with the property that 3?(V(x)) is uniformly bounded from below on all the chosen contours and tends 
to oo polynomially (in the length parameter) on the contours that extend to z = oo. In more detailed 
terms: 

(a) Consider a pole z = c of V of order k > 2: we attach to it k — 1 "petals" approaching z = c along 
asymptotic directions in the sectors where $t(V(x)) — > +oo. We also attach a "stem" extending to 
oo and asymptotic to a direction such that $l(V(x)) — > oo. 

(b) For a simple pole z = c of V' , if the residue is a positive integer (i.e. e~ v has a pole at z = c) we 
choose a small loop around the point, if the residue is a negative integer we take a contour from 
z = c to oo, if the residue is non integer we take a loop coming from oo and returning to oo (with 
the same restriction as above for the asymptotic direction). 

(c) We choose also arbitrary segments joining a certain number of points z = a to oo (along admissible 
directions). These latter contours are called "hard-edge" contours because the pseudo measure 
d/i = e~ v ( z 'dz has a limit at z = a and integration by parts yields a boundary term. 




Figure: The contours for a typical semiclassical moment functional. Here V'(x) 
has a pole of order 4 at oo, of order 4 at C3 and simple poles at C2, C3 with nonin- 
teger and negative-integer residue respectively. The contours originating from 
the di's are "hard-edge" contours. The shaded sectors represent the asymptotic 
"forbidden" directions for approaching a singularity. One of these sectors at 00 
in Figure does not have a contour surrounding it because such a contour would 
be "homologically" equivalent to minus the sum of all others. 
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7.1 Differential equations 

We first analyze in this situation the infinite-dimensional differential equation that the BOPs satisfy. 
The natural differential operation in this setting is not d x but rather xd x . Using the recurrence relations 
involving multiplication by x and the orthogonality relations 



a;p(a;) = Qp(x) ; xr t (x) = r t (x)Q 
xp'(x) = Dp(x) ; xr'(x) = Dr(x) 



pr'e v dz 



we can obtain the following identity by integrating d z (zpr t e y W) 



D 



+ £>*- (zp(z)r t (z)e" v(z) ) 



QV'{Q) 



(7-3) 
(7-4) 

(7-5) 



(7-6) 



The notation that we now adopt is that J stands for the linear combination with coefficients Xj of 
integrals on the oriented contours Tj and the evaluations stand for the evaluations at all endpoints 
of the given contours, multiplied by the corresponding coefficient k and the appropriate sign according 
to the orientation. The matrices D and D are lower triangular and on the main diagonal they can be 
explicitly computed 



xp n = np n + previous ; 

xr' n = ( (— n)(l - £n+i) + i n +i ) r n + previous 



x— r* n = nl n +ir* + previous 



This implies the following Virasoro scaling constraint 

{QV'(Q)) nn + (z Pn r n e- v ) 

Note that we also have 
n-l 



= 1+4 



+i 



n£ n+1 



k=0 



J2 (Q v '(Q))kk + {zp 



kr k e 



f 1 + + - n{£ n + 1 



(7-7) 
(7-8) 

(7-9) 



(7-10) 



(7-11) 



k=0 



The parts of D, D below the main diagonal are now expressed in terms of Q and the boundary terms 
only 



D< = (QV'(Q)) < + (z(p(z)r t (z)) < e- 1/ ( z ) 
D< = (Q t V'{Q t )) < + (z(r(z)p t (z)) < e- ,/ ( z ) 



(7-12) 
(7-13) 



Note that -below the main diagonal- the matrices D and D are of the same form as the deformations we 
were considering previously; more precisely they correspond to a variation by F(z) = zV'(z) and a linear 
combination of variations supported at the endpoints of the contours Tj . The folded version of this ODE 
can be obtained from the formulae (|6-91 16-141 16-211 16-23(1 with the only modification that comes from the 
diagonal part of D. Using I|7-1U|) for the diagonal part the reader can check that the result is 



n 
xV'{x) - 1 



In 



-W, 



■W n ,n-l 
(n-l)*,n-l 



W n , n * 



zc 



-V(z 



In 



-p n r n -i 



• n-V 



PnPn 
r n -lPn 
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W := VqxV'(x) = 



QV'(Q)-xV'(x) 



We remark that the last "boundary" term consists of simple poles with nilpotent residues. 
For the full matrix \ n the differential equation is 

xd x Xn(x) = V n {x)xn{x) + Xn(x)T> R (x) 



(7-14) 



(7-15) 



V R (x) 



xV'(x) - zV'(z) c -y( z ) dz . 



-V(z) 



-xV'(x) 



(7-16) 



Together with the differential equation and the deformation equations we recall that we also have differ- 
ence equations 



Xn 



R n (x)Xn-l 

x-f3„ 



n > 1 



R n (x) = < 



(_ljn+l 

x-/3 n 



(-1) 



Tt + 1 



Ur, 



if t. 



if in = 



(7-17) 
(7-18) 



The ladder matrices R n are simply obtained from the transfer matrices H3-2|) by using the normalization 
of the polynomials as in (|4-5H . We have thus proved 

Theorem 7.1 The matrix Xn satisfies the following system of difference- deformation- differential (DDD 
for short) equations 



Xn = R„(x)Xn- 



da; 



(7-19) 
(7-20) 

Sf Xn = U^> X n+XnU^>' K (7-21) 

where f denotes either any function or formal power series provided that C(f(z)z k ) is well defined for 
k G Z or any derivative of the Dirac delta function supported at any point a^O. 

We observe that the right action of the differential-deformation equation is independent of n. This 
suggests that we can perform a "right gauge" change to dispose of this part. Indeed we define the new 
object r„ which will be the focus in the rest of the paper 



1 







3 V(x) 



-V(z) 



-dz 



(7-22) 



It is easy to verify that this change of gauge eliminates the right-actions for the differential equation and 
for any deformation of V(x) and/or the endpoints of integration. The first column of T n is the same 
as the first column of Xn and hence contains the LOPs. The second column contains now the following 
auxiliary functions 



-dz 



-V(z) 



x n-i-e ne v(*) f r n -i(z)e- ^ 



(7-23) 
(7-24) 
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We note that the auxiliary functions are piecewise analytic functions off the contours Tj: it is a matter 
of routine inspection to read-off the relevant Riemann-Hilbert data. We defer this inspection to a later 
section. 

In terms of the matrices r„ we have a DDD system of more standard form, without right multipliers. 



Theorem 7.2 The following system of Difference- Differential- Deformation equations is Frobenius com- 
patible 

r n = R n {x)T n _ x (7-25) 



x— r„ = V n T n (7-26) 
da; 



_d_ 

da; 

5 f T n = U<pT n (7-27) 

where f is as in Thm. j7.1\/ . 

A few remark are in order here: by choosing / in Thm. 17.2fl to be an ordinary function one can vary 
the potential V by V — ► V + ef and hence all flows of the generalized Toda hierarchy are here included. 
However we can also choose / as a distribution or linear combinations thereof. Clearly if we choose 
the point a arbitrarily outside of the singularities of V{x) we still have a compatibility of the resulting 
system but we will change the structure of the singularities of T) ni which falls outside of the standard 
theory of isomonodromic deformations. For example, adding a S a corresponds to adding a term ln(x — a) 
in the potential and adjoining a small circle around a to the set of contours Fj's. 

Vice- versa the cases in which / is a distribution which does not alter the singularity structure of T> n are: 

1. Movement of the endpoints which contribute to the boundary term 6 : then we have 

/ = ±xe- v{a) S a (7-28) 

where the coefficient k is the coefficient of the contour Tj which has a as endpoint and the sign 
depends on the orientation of Tj . 

2. Movements of poles of order k (if any) of the pseudo-measure e~ v dz: then we have 

f{a) = ±^k5^ +1) {z)e- v ^ z) (7-29) 

where the coefficient x is the coefficient of the loop encircling a, the sign is chosen according to the 
orientation of the contour and V r {z) is the part of V which is regular at z = a. 

8 Spectral curve and Isomonodromic tau function 

The objective of this section is that of expressing the spectral curve of the connection d x — -T> n (x) 
in terms of the logarithmic derivatives of the Toplitz determinants; this will be the essential bridge to 
connect with the isomonodromic tau function in the coming sections. We prove the following theorem 

Theorem 8.1 The following formula holds 

det (yl - \v n {x)j = y 2 - (v'(x) + ^f)v + 

+ I T r„ (QXM^XM.) + I f^VlWN 
x \ Q — x ) x V x — z ) 

L n := n - 1 - 4 , (8-2) 



They corresponds to those endpoints of the contours Tj for which . e 7^ 0. 
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where II„_i = diag(l, 1, . . . , 1, 0, 



(n nonzero entries). 



Before proceeding to the proof we remark that this formula would be valid for an arbitrary smooth 
potential; quite clearly, however, in this case the spectral curve would not be an algebraic curve. 
Proof. 

We need to compute the two spectral invariants of the connection; the main tool is to use the compatibility 
between the ladder relations and the connections T> n (x). Indeed from the compatibility between the 
difference and differential equation and from the explicit expression for T> n (x) (|7-14|l we can express 
recurrence relation for the spectral invariants of T> n (x) . The trace is computed by sight 



Tr(p n (x)) = xV'(x)+n-l-e n 
From the compatibility of difference-differential equations we have the gauge property 

T) ri —\ Rn T-^nRn Rn R71 

The gauge term is explicitly computed to be 



Rn T^ii Rn 



X Rrr Rn 



T) n -\ + & Rn l R' n 



i-iy 



7n-l 
1 





1 






if tr. 



if tn = 



(8-3) 



(8-4) 



(8-5) 



(8-6) 



These formulae imply a recurrence relation for the quadratic invariant. 

Tr(££)) = Tr^ O + 2 Tr (V^xR^ 1 R' n ) + T^xR^ 1 R' n ) 2 ) 



-7) 



For the line case i.e. £ n = 1 and using the form of the recursion matrices R n together with the fact 
that in this case r*_ 1 = (— l) n ~ 1 p n ^i we find 



Tr(Vl) = Tr(^_ 1 )_ 2a 



/ QV'(Q) - xV'(x) 



Q - x 



2x 



-V(z)p„_ir„_ 



x — z 



For the circle case £„ — instead we have 



Tr(V 2 n ) 



Tr(^_ x ) + 2 (n - 1) + 2 7 „_i 



-W n _ i.Ti-2 H ; = W n -i t (n-iy 



-27„_i 



-V(z) 



Pn-l 



-r n -2 + 




(8-9) 



Using the identity (|2-7|l together with the definitions of the biorthogonal polynomials and the various 
normalization factors l|4-5|) one can see that 



r„-2 + 



=Pn-l = 



7n-l 



-r n -x 



(8-10) 
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and hence 



(-1)"A £ "_ _1 1 
W„_i,„_ 2 H : — p==W„_i i( „_i)* = : £ z | zlfp„-ir n _i 



Therefore the recursion for the circle case is 
Tr(^) _ TriVl^) = 2(n - 1) + 1 - 2 



7n-l 



Q(Qn<2)-^'(z)) 



(8-11) 



Q - a; 



- 2 

n — l.n— 1 



z 2 e y(z) p„-ir„_i 

a; — z 



= 2(n - 1) + 1 - 2 (q^'(Q)„-i,„-i + C*e- v <*W-irw)| e3< ) + 2xV'(x) + 

V Q-x /„_i,„_i 

= 2(n - 1) - 1 - 2£„ + 2a;V (a;) - 2x 



Ox 



zc y(z) p„-ir„_i 



a; — z 
/ QV'(Q)-xV'(x) 



71—1,71—1 



-2a; 



ZC ^^^n-lT-n-l 

a; — z 



(8-12) 



9;* 



Summarizing, in the two cases we have found 

Tr(Vl) - TriVl^) = 2 (xV'(x) - l n + (n - 1) - t ) ( | . . /, . , + 



2x ( QV(Q)-xV'(x) 
V Q-x 



- 2x 

71—1,71— 1 



zc v{z) y n -\r n - 1 
x — z 



(8-13) 



To complete the computation we need to find Tr(2?i ) or -equivalently- det(2?i). We have 

-l 



dct ( -Di ) = det 

. x 



P'l V4 



r* ® 



= det 



Jne~ v(a ° det 

V'(x)C z f — 

\x - z 

y'(Q)-y'(aO 

Q - a; 







This implies 

detX»i(a;) = x 2 



v'(Q) - v'(x) 
Q-x 



oo 



oo 



-V(z) 



08' . 




Z^'(z) 




— z 









P-V(z) 



-V(x) _ 



Pprp \ 



(8-14) 



Por \ 



x — z J 



=0 



= x (V(Q) 00 + (poroc-^))!^) + s ( gf^ 
(QV'{Q) -xV'(x) 



V Q - a; 



oo 



V Q - a; 

a; — z / 



oo 



ze ^^Vo^o 

x — z 



(8-15) 



9;* 
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Hence (4 = 0) 



TrCDi 2 ) = (xV'(x)f - 2x 



( QV'{Q)-xV'{x) 
V Q~x 



2j { ze-n*) poro \ 
x — z J 



(8-16) 



Tr(V n 2 ) = (xV'(x)) 2 + 2xV'(x) (n - 1 - i n ) - (n - 1 - 4) + 2^(fc - 1 - 4)(1 - 4) + 



fe=i 



Q - x 



x — z 



(8-17) 



Using this expression for the quadratic invariant we can obtain the following formula for the characteristic 
polynomial 

det(p - X> n (a;)) = y 2 - (xV'(x) + n - 1 - £ n )y + K n + 

+xTin (QV'(Q)- X V>( X )\ + x /ze-vWp^-xr 
V Q — a; / V x — z 



(n-l-e n )(n-e n ) 



^(4 + i-fc)(i-4) 



(8-18) 



fc=2 



The last crucial observation is that K n = for all n: this is non-obvious at first sight and it is true only 
because £ n is a weakly increasing sequence of integers. Indeed one can check that 



Kn+i - K n = -4+1 (1 - 4+i) , 



(8-19) 



so that Kn+i = K n = K\ = 0. To conclude the proof we note that the spectral curve of Ij8-18l) is simply 
related to that of the connection by y = xy. This ends our proof. 



9 Isomonodromic deformations 

By Thm. 17 .21 we have compatible systems of Difference, Differential, Deformation equations 



r„ 

d x T n 
S f T n 



—V r 

X 



(9-1) 
(9-2) 

(9-3) 



The compatibility of this system entails isomonodromic deformations for the connection d x — ^.D n . Note 
that this connection has the same singularity structure of V'(x). In order to have isomonodromic defor- 
mations in the sense of Miwa-Jimbo-Ueno we need to impose that V'(x) be a rational function. Then 
the deformations of V(x) which give rise to the setting in MJU are those which do not alter the singu- 
larity structure of V(x); this is why the most general setting compatible with this requirement is that of 
semiclassical moment Junctionals. 



9.1 Spectral residue-formulae 

Mimicking the approach of ^ we can express the logarithmic derivatives of the shifted Toplitz deter- 
minants A^™ in terms of residue formulae involving the differential ydx on the spectral curve defined in 
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Ln, \ 1 If v , L r 



Thm. O 

Y±(x) := ~ ( + ^ ) ± i J ( + ^ ) " (9-4) 

(9-5) 



1 ^V"(Q)-a;rW\ , 1 fze-^p'IIjv-ir 
Pix) = -Tr„ 



V Q — x 



Indeed we have 



Theorem 9.1 Let V'(x) be rational. 

(i) Suppose that x = c is a pole of order d + 1 



(c) 



V W = E T7TZ37 - 4 C) - c) + O(l) (9-6) 

,7=1 ' 



t (c) 



.7=0 V ' 7 



TTien we /lave 



i ( 7 c) = - res F + (a;)(x - c) J dir , J = 0, . . .d (9-8) 

x—c 

d In A^" 1 

-iesY^{x){x - c)' J dx , J=l,...,d (9-9) 



at W J x = 



to 



C?C z=c \ ^— ^ (x — c) 

\J=0 X ' 

(ii) Suppose that x — oo is a pote o/ V with degree d, namely 



d+l .(oo) 

y(ar) = E^r 2;J + °( ln2: ) ( 9 " n ) 



d+1 



y'(a;) = ^^ oo) x J - 1 +0(l/x) (9-12) 



,7=1 



TTien we /iaue 



= -^rc^r + (.x)a;" J da; , J=l,...rf+1 (9-13) 
<91nA £ ™ 1 

» = - res x J Y_(x)dx , J=l,...,d+1 (9-14) 

()t ^ £=00 

(iii) Lei x = a be a hard-edge 7 , namely a point of the boundary of one of the contours {Tj} such that 
\V(a)\ < oo. Then 

= 1 res ±Tr(D n fdx (9-15) 

oa 2 x=a x z 



7 This means that this is one of the points contributing to the boundary terms. 
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(iv) Finally we have 



res Y+(x)dx = L n = n — 1 — l n — J]] to 

c= finite pole of V 



(c) 



res Y + (x)dx = £ n + 1 + 

x—oo 

Proof. We start by noticing that 



-16) 



0{{x-c) d+1 ) for case (i) 



Y 



V'{x) + 



V(x) 



V'(x) + ±* T-r~Y + 0(x~ d - 2 ) for case (ii) 



(9-17) 



At this point formulas I|H-8I9-13I l9~THjl follow immediately by noticing that V/(V'(x) + L n /x) = 0(1) in 
all cases and by straightforward computation of residues 8 . As for the remaining formulas we have, for 
case (i) 



res (x — c) J Y- (x)dx = res (x — c) 



(x — c) J 
res : r ., ; ' - Tr» 



V'(x) + ^f x=cxV'(x)+L n 



xV\x) - QV'(Q) 



n-l 
k=0 



(x-c)- J xV\x)-zV'{z) . 

res T777 n T t Pn(z)r n (z) 

x=c XV (X) + L n x — z 



-Tr„(Q - c)" J = Jfl w ln A<" , J = 1, .. . ,d , 

j 



(9-19) 



and similar computation for the c-derivative. Here we have used the formulae (|5-12|l expressing the 
variation of ln A^™ under an infinitesimal deformation of the type ensuing from an infinitesimal change 
of the parameters tj. 

For case (ii) the computation is completely parallel except for the last J = d + 1 residue. Indeed 

fxV'{x)-QV'{Qy 



V(x) 

res x J Y_(x)dx = res x J T 



res 



fc=0 



res 



00 V (x) + *f x =°° xV'(x) + L, 
x J xV'Cx) - zV'(z) . . , , 

-Pn(^)r n (2J 



-Tr„ 



ie=oo + L n X—Z 

= -Tr n Q J = Jdjoo) ln A*" , J = 1, . . . , d 



(9-20) 



For J = e?+ 1 one has to use a similar manipulation but has to use the refined asymptotics <|9-17[) : indeed 
we have 

res xV'{x)Y^{x)dx = -Tr n QV'(Q) + ( n 2 - nL n - Tr n QV'(Q) - 0e" v(z) p*n n _ 1 r 

= -Tr n QV'(Q) (9-21) 

+ GO" 2 ) . (9-18) 



8 Note that at infinity 



V(x) _ n 
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where we have used l|7-lll) together with the definition of L n = n — 1 — l n . This proves, together with 
the residues (|9-20j) 

res x d+1 Y^{x)dx = -Tv n Q d+1 = (d+ ln A*° . (9-22) 

Finally, for the case (iii) the computation is immediate using the formula for TrU^ I|8-17f) . Q. E. D. 

10 Riemann— Hilbert problem, Tau function 

Direct inspection of the asymptotic behavior of the biorthogonal polynomials and second kind functions 
allows us to ascertain the Riemann-Hilbert data for this problem. We start by noticing the following 
formal asymptotic behavior of the auxiliary functions entering in r„ 



-V{z) 



Pn(z) 



x — z 



= < 



(_)» s -*»-a e v(*y _j±!(i + 0{x- x )) for x -> oo 

_ x n-l-e ne V(x) A n+1 n+HM) forz^O (10_1) 



near poles of V'(x) 



IJ n-l 



i-<„ p v(x) f e v ^r n _ 1 (z) 



x — z 



= < 



x -t»-i e vw 






x n-l-in e V(x) 


V 1 (1 + 










e v( - x] ^/%c n - ( 


™ + (Q ~ c )o,n-l 





where we have used the definition of the LOPs (|4-5|l and the facts that 



r n _i cx z 

This implies the following formal asymptotic data for T n near all the singularities. 
At x = we have 

r„(aO ~ Gi°) 



1 

x n ~ 1 ~ e ' l e Vairi3 ' ^ 



(10-2) 

for x —>■ 
near poles of V'(a;) 



(10-3) 
(10-4) 



(l + 0(x)) 



(10-5) 



At x — oo we have 



'A "A 



r^z) = g,[ ! oo) 



detGW = 



a;" 

a; - f n- 1 e v 'sin 9 ,oo(z) 



1 + 



(10-6) 



(10-7) 
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Near any other pole x — c of V '(x) we have 

r„(x) = Gi c 



Q(°°) 



1 o 

easing, e(x) 



y/hn-1. 



(H-O(x-c)) 



(10-8) 



(10-9) 



G (c) 

v * r» 



Pn(c) 



l ^e g , C ( C ) 



(10-10) 



where in all these formulae the notation V S i ng _ p (V re g tP ) denote the singular (regular) part of V at the 
point p. 

Near a hard-edge point i = «we have 0] 



G (o) 



1 ±xln(a; — a) 
1 



(1 + 0(1 -a)) 



G (a) := 



Pn (a) e^W£( p " ( ^" (a) ) 

r* , (a) a n-l-^ e V(a)£ / r„(s)-r„(,) 
ro— 1 V / \ a— z 



(10-11) 



(10-12) 



Together with these data we also have the jumps across the contours 1^ defining our moment functional: 
the situation in this respect is identical to 0]. In essence the matrix r„(x) has the following jumps across 
the contour Tj 

t, f \ y-i / \ [ 1 2z7T>f,; 

T n (x) + = T n (x)- 







1 



(10-13) 



Note that these jumps can be interpreted -depending on the point of view- as the Stokes' multipliers of 
the problem near the singularities. 



10.1 Isomonodromic Tau Function 

Using the results of 0j we find that the Jimbo-Miwa-Ueno isomonodromic tau function is given by the 
same differential formulas in Thm. l7.2l provided that we substitute the spectral curve of the connection 
with the spectral curve of the connection in the traceless gauge. In our situation the trace of -T) n {x) is 
V'(x) + n ~ 1 ~ e " n so that we perform a scalar gauge transformation 



A^ MU) = l -V n (x) - \ (v'(x) + 



n-l-L 



L2x2 



(10-14) 



This implies that the spectral parameter ujmu h &s the following relation to the previously employed y; 

(10-15) 



Vjmu = V + g I v K x ) H 

This in turn implies that -up to multiplicative factors independent of the isomonodromic times- 



F[V]tjmu 
\nT[V] = -~ 



V'(x) := V'{x) + 



Y] res (x)V reg , c (x) 

£ — ' x— c 3 ' 

finite pole of V' 

n - 1 - i n 



(10-16) 
(10-17) 

(10-18) 
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where Vg ing c ( V re g,c ) denotes the singular (regular) part of V' at the pole c. 

Example 10.1 For example let us consider the case relevant to the problem of the probability of the 
longest increasing sequence of random letter in a word of fixed length 1161/ 

M ( - \ 

V(x) = -tx-J2 ka ln ( ) (10-19) 

In this case a direct computation (£ n = 0) gives for T the following expression 

M M M M , s kck/3 

^ = -\ E k ^ + V 1 E ka H ~ rJ) + 1 E k " ln (- r « 2 ) - \ E E ln rj y^ ( 10 - 2 °) 

which is the result obtained also in formula (3.76) : note that in that formula ro = and ko = 
n — ^] . r Q and a ZittZe of algebraic manipulation shows the equivalence. Moreover the signs inside 
the logarithms in \1U-2U\) are in fact irrelevant since omitting them would amount to multiplying T by 
a constant independent of the isomonodromic times, and hence could be reabsorbed in the definition of 

TJMU ■ 



11 Schlesinger Transformations 

From the asymptotics that the shift implemented by the matrices R n are -in the language 

of isomonodromic deformations- what is known as elementary Schlesinger transformations. Specifically 
the shift n i— > n + 1 corresponds to the following two types of elementary Schlesinger transformations 
according of the type of move (circle or line) (refer to formulae 110-61 and I10-8JI : 

Circle move. The Schlesinger transformation adds one to the first entry of the formal monodromy 
at oo and subtracts one from the second entry of the formal monodromy at zero 

Line move. The Schlesinger transformation adds one and subtracts one to the first and second 
entry (respectively) of the formal monodromy at infinity, leaving the formal monodromy at zero 
unchanged. 

However we can obtain a third type of elementary Schlesinger transformation by considering two distinct 
sequences of LOPs corresponding to two (weakly increasing) sequences of {4i}'s. Suppose indeed that 
we consider another sequence of LOPs and the ensuing connection xd x — T> n (x) for some fixed n where 
the only difference between the two pairs of LOPs is that one (or more) circle-moves have been replaced 
by a line-move (or vice-versa) along the chain for n' < n: the only difference in the formulas will be that 
t n — (-n i 1. This is implemented by the "circle-to-line" transformation T„ 12-151 fsuitablv normalized). 



Pn 



b c 
a+- - 

x x 



d 

x 



e 
x 



Pn 



(11-1) 



where the coefficients a,b,c,d,e above can be obtained explicitly in terms of shifted Toplitz determinants 
using the form of T n l|2-15f) and the normalizations 1)4-5(1 . and the polynomials p n ,^n-i refer to the 
elements of the sequence of biorthogonal polynomials associated to the sequence {£k}- such sequence 
differs from {£k} because l n = £ n — 1, namely there is & k$ < n such that 1^ = — 1, Vfc : ko < k < n. 
We therefore add the following third type of transformations; 
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Circle-to-line move. The Schlesinger transformation subtracts one to the second entry of the 
formal monodromy at oo and adds one to the second entry of the formal monodromy at zero. 

This last type of transformation shows that the orthogonal polynomials on the line and the orthogonal 
polynomials on the circle are related by a sequence (n — 1) Schlesinger transformations and at each step 
the Laurent biorthogonal polynomials that are obtained are those appearing in the solution of integrable 
lattice hierarchies associated to elementary orbits 9 . 



12 Conclusion 

As a general "philosophy" , it is acknowledged in the literature that KP tau functions and isomonodromic 
tau functions are often, if not always, related one to the other, in the sense that a KP (or Toda) tau 
function is an isomonodromic tau function for a suitably chosen isomonodromic deformation. In the case 
of orthogonal polynomials this relation was explored in JHj for some class and extended in [HIE]- In this 
paper, this relation has been confirmed once more for the particular generalized Toda systems associated 
to "nonstandard" minimal orbits of the Borel subgroup: the natural bridge between the Hamiltonian 
and isomonodromic treatment is provided by the solution of the inverse spectral problem in terms of bi- 
orthogonal Laurent polynomials. It is to be expected that, whenever a description or formulation of an 
integrable dynamical problem in terms of (bi/multiple-orthogonal) polyomials is available, then a suitable 
definition of the tau function for the associated isomonodromic problem should tie the Hamiltonian tau 
function with the isomonodromic one. For instance, in the case of the biorthogonal polynomials arising in 
the study of two-matrix models [5][7] a natural isomonodromic deformation of a polynomial connection can 
be derived; however the connection is a highly resonant one and at present a definition of isomonodromic 
tau function for resonant deformations of connections is not available. However it is possible to formulate 
such a notion [H] and the connection can thus be positively be established. 

As it is recalled in the appendix to follow, the Laurent orthogonal polynomials which we have investi- 
gated in the present paper are related to the solution of the inverse spectral problem for Toda-like systems 
associated to certain minimal (or elementary) irreducible orbits. There exist in fact other minimal orbits 
for which a treatment in terms of orthogonal polynomials of some sort is not readily and generally avail- 
able, although inspection of specific examples leads to expect that it is possible to overcome the difficulty. 
It is our intention to pursue the topic in future publications. 



13 Appendix: minimal irreducible co-adjoint orbits 

As it was mentioned earlier, every n x n principal submatrix of the Hessenberg matrix Q that defines 
recurrence relations (5-1) belongs to a 2n — 2-dimensional co-adjoint orbit of the Borel subgroup B„ of 
invertible upper triangular matrices in sl(n). However, not every low-dimensional co-adjoint orbit can 
be obtained this way. In this appendix, we give a description of all irreducible co-adjoint orbits of B„ in 
sl(n) that have a minimal dimension In — 2. 

First, we introduce some notations. Let b_ be a subalgebra of lower triangular matrices in sl(n). 
Denote by J an n x n shift matrix (Is on the first superdiagonal and Os everywhere else) and let Hess n = 
J + b_ denote a set of lower Hessenberg matrices. An element Q S Hess n is called reducible if it has a 

Qn Q12 



Q 2 



block upper triangular form Q = 

irreducible otherwise. 

Orbits of the co-adjoint action of B„ on Hess n are given by 



, where Qn is a k x k matrix ( < k < n ). Q is called 



Oq q ={J+ (Ad b Q )<o : b e B„} . (13-1) 
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It is easy to see that if Oq contains a reducible (resp. irreducible) element then every element of Oq 
is reducible (resp. irreducible). Therefore it makes sense to talk about irreducible orbits of the co-adjoint 
action. Our main goal in this appendix is to prove the following 

Theorem 13.1 An irreducible co-adjoint orbit o/B„ in Hess n has a minimal dimension (2n— 2) if and 
only if it contains an element Qq of the form 



Q = J + H + J2 E . 



i a ,i a -i— e a -i 



a=l 



where 

1. Si £ {0, 1} and £o = 

2. 1 = io < h — e\ < it < h — £2 < «2 < 

3. H = 2_j h a E aa . 

a£{l,...,n}\{i ,...,ik-i} 

1 

: * 



< h-i - £fc-i < ik-i < h 



Qo 



1 

1 

★ 1 

: * 



1 



1 



1 



1 

1 



1 



(13-2) 



(13-3) 



[An example with n = 11, k — 3, i\ = 4,i2 = 8,13 = 11, e\ = 63 = 0, £2 = 1./ 

Remark 13.1 The case of H = hE nn and e a — 0(a = 1, . . . , k — 1) was studied in 1 10f . It is orbits 
of this type that can be studied via associated LOPs of the type appearing in this paper. Note that, in this 
case, parameters £j that were used in the main body of the paper are related to i a via 

£j = max{i Q : i a < j} . 

An investigation of propeties of moment functionals connected with a more general minimal orbits de- 
scribed in Theorem \l'S.l\ will appear elsewhere. 



Define a staircase pattern (I, e) as a collection of pairs of indices 

(I, e) = {(it, 1), (i 2 , i\ - £1), . . . , (ifc = n, ik-i - £fc-i)} , 

where 

1 = i < ii - £1 < h < i 2 - £2 < «2 < ■ • ■ < ik-i - £fc-i < ik-i <ik = n . 



(13-4) 
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In what follows we will often use a notation 

Ja ^a — 1 ^a— 1 ■ 

We say that Q £ Hess n has a staircase pattern (/, e) if 

Qi a ,j a 7^ and Qy =0 for i > i a , j < j a+ i (a = l,...,k) . 

The set of all matrices in Hess n that have a staircase pattern (I, e) will be denoted by Hess(I, e). For 
example, if I = {2, 3, . . . , n} and e = {0, 0, . . . , 0}, then Hess(I, e) coincides with the set of n x n Jacobi 
matrices. An immediate property of the set Hess(I,e) is that it is stable under the co-adjoint action of 
B„, since corner entries Qi a ,j a and the entries "under the staircase" Qij — i > i a , j < j a +i have only 
zeroes to the left and below and, thus the former are being acted upon only by the diagonal part of B r , 
and the latter cannot be made non-zero by the co-adjoint action. 

Let us fix a staircase pattern (I, e). To begin the proof of Thcorcm ll3.ll we first employ the strategy 
used in to study generic staircase orbits. 

Lemma 13.1 If Q G Hess(I, e) then there exist Q € Oq such that 

Qiaja 1 

Qija = (ja < i < ia) 

Qi a j = (ja < j < ia and j ^ jp : (3 < a , jp < i a ) 

(a = l,...,k). (13-5) 



Proof. First, we use a diagonal conjugation to reduce Q to an element with all corner entries equal 
to 1 : Q — > Ad* D Q = D^ 1 (Q — J)D + J, where D = diag(di, . . . , d n ) = Dk ■ ■ ■ D\ with diagonal matrices 
D a defined by 

(n \ -/ 1 { ^ la 

\u a)ii - | d ^ = ( D -i^ D -i QDl . . . D a -i) iaja i = i a 

Next, we use the co- adjoint action induced by a sequence of elementary upper-triangular matrices (each 
depending only on one parameter only) to set as many as possible of the entries in rows and columns 
occupied by corner entries equal to zero. More precisely, to eliminate an (z,j a )-entry (J a < i < i a ) 
using the corner entry (i a ,j a ), one employs ^4c?^ 1+ g £ y Similarly, to eliminate an (i a , j)-entry 
(ja < j < ia), one uses Ad*^ 1 _g, , E . ,y Note that when we write Qij a (resp. Qi a j), we refer to entries 
of the "current" value of Q, i.e. to the element that belongs to the orbit through the initial Q and that 
has been obtain through the sequence of transformations already applied. 

The order in which we apply these elementary transformations is defined as follows: we first set 
to zero the entries in the 1st column (going down the column), then in the i\st row (moving right), 
then in the j2nd column (moving down), then in the ^nd row (moving right) etc. Through the entire 
process, we want, for every I < to, to use an elementary matrix of the from 1 + xE[ m at most once. 
This means, in particular, that any (i a , t 7^)-entry, where < a and jp < i a cannot be touched, since a 
matrix of the form 1 + xEj fj i a has already been used to eliminate the (jp, j a )-entry. This explains why 
entries {(i a ,jp) : (3 < a, jp < i a } are excluded from the list of entries in l|13-5|) . On the other hand, all 
non-corner entries that are in the list can be set to 0, regardless of their initial values. Q.E.D. 

Corollary 13.1 For each Q € Hess(I , e) the matrix entries specified in \ 13- 3(1 are independent functions 
on Oq. 

Proof. It suffices to notice that applying to Q constructed in Lemma ll3.1l elcmcntarv transformations 
of the same type that was used in its construction, but in the reverse order and with arbitrary parameters, 
one can obtain an element in Oq with arbitrary nonzero values of the corner entries and arbitrary values 
of non-corner values specified in Ijl3-5|l . Q.E.D. 
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Lemma 13.2 If, for some 1 < a < k, e a > 1, then, for any Q £ Hess(I,e), dim0Q > 2n — 2. 

Proof. Denote by M(I, e) the set of pairs of indices that appear in the list given in l|13-5ll . In view of 
the corollary above, we only need to show that, under conditions of the Lemma, the number of elements 
in M(I,e) is greater than 2n — 2. We will also show that, if < e a < 1 for a = 1, . . . ,k — 1, then 
#M(Z,e) = 2n-2. 

We will use an induction on k and n. Clearly, if k = 1 then £o = and #M (/, e) = 2n — 2. Moreover, 
#M(7, e) = 2n — 2 for any fc, provided e Q = for all a. Let now k — 2 and ei > 0. We are looking for a 
number of elements in the set {(1, 1), . . . , (it, 1), (it, 2), . . . , (it, it — 1); (i\ — £i, H — £i), . . . , (ii — — 
ei,«i — ei), (it + 1, U — £i, U — £i)> • • • j ii — £i)> • • • (n, it — 1), (n, it + 1), . . . , (n, n — 1)}, which is equal 
to2(i 1 -l)+2(n-(ii-ei))-2 = 2(n + ei -2) | =^»-2 if ei = 1 

For k > 2, let s be such that j s < ij < j s +i- We first consider the case when there is no r such that 
j r = it- Then the set M(I, e) \ {(1, 1), . . . , (h, 1), (it, 2), . . . , (it, it — 1)} has the same cardinality as a set 
M(I',s'), where 

i 1 '^') = {(*2 ~h> 1)) (*3 - J2,J3 - i2 + 1), ■ ■ • , -J2,is - J2 + 1), (*a+l ~ J2, ~ J2), ■■■,(«*; - J2, jfe 

{(«2 - il + £l, 1), (13 - h + El, «2 - il + £1 - (£2 - 1)), ■ ■ ■ , (is -h + £l, - il + £1 - (e s -i - 1)), 
(i a+ i - ii + et,is -ii+ei- £ s ), . .., (i k - it + £i,4-i -«i +ei - £fe-i)}, that is 

n' = 4-i = «fc - h + £1 , 

i' a = i a +i - it + £1 a = 1, . . . , A; — 1 

and 

£^ = £ Q+1 - 1 (1 < a < s - 2) , e' a = e a+1 (s - 1 < a < k - 2) . 

If s > 2 then £i > 1 and £ s = i s — j s > i s — it > s — 1 > 1, so that s' s _ 1 — e s > 2 and, by the induction 
hypothesis, #M(/',£ / ) > 2(n - ii + ei - 1) > 2(n-ii) and #M(7,e) > 2(«i - 1) + 2(n- ii) = 2(n-l). 
If s = 1 then £^ = £ Q +i for 1 < a < k — 2 and 



#M(l',e') 



= 2(n-ii+£i-l) if all e' q < 1 
> 2(n — «i + st — 1) if some e' a > 1 



and thus, #M(7, e) = 2(ii — 1) + #M(/', £ 7 ) is greater than 2n — 2 if e a > 1 for some a > 1 and is equal 
to 2n — 2 otherwise. 

Finally, consider the case when j r = i\ for some r > 1. If r > 2 then £ r = i r — j r = i r — i\ > r — 1 > 2. 
Define (I,e) = (I,e) \{(t2,j2), ■ ■ ■ , (V-i>ir-i)}- Then (/,£*) still defines an irreducible staircase pattern, 
fe = #(/, e), fc and #M(I, e) > #M (J, e) > 2n - 2 by the induction hypothesis. 

If r = 2 then e x = 0, j 2 = i x and #M(J,e) = 2(i x - 1) + #M (/',£'), where (/',£') = {(i 2 - i x + 
1, 1), (?3 — ii + 1, %2 — it + 1 — £2), • • • j (n — it + 1, (ifc-i — «i + 1 — £fc-i)} and, again by induction, the 
statement follows. Q.E.D. 

We are now ready to complete the proof of Thm. 113.11 

Proof of Theorem 113.11 Assume that dim0Q = 2n — 2. We have shown that if Q G Hess(I,e) 
then e a < 1 for a = 1, . . . , k — 1. Assume that the latter condition is satisfied and consider the element 
Q constructed in Lemma |13. II Suppose that some non-corner entry Qij (i > j) is nonzero. Then, by 
construction of Q, j 7^ j a (a = 1, . . . , k). Define a diagonal matrix D by 



D u 



1 if I ^ j or I ^ jp or i ^ ^ 

g? if I = j or ( I = jp and i = ip ) 



Then Ad^Q = D l QD has the same values as Q in the entries specified by (|13-5|) but (Ad* D Q)ij = 
d~ x Qij- This means that the matrix entry Qij viewed as a function on Oq is independent of the matrix 
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entries specified by Ijl3-5|l . which is in contradiction with dim0Q = 2n — 2. Therefore, Qtj = for all 
ihj) 7^ (*q, j^). Since, by Lemma |13. II Qj a j a — for a = 1, we proved that dim0Q = 2n — 2 

implies that Oq contains an element of the form l|13-2|) . 

To prove the converse consider an element Q defined by (|13-2|) . Clearly, for any b 6 B„, Ad^Qo = 
Ad^(Qo — if) + H, therefore it is sufficient to consider the case where h = 0. In other words, we are 
interested in parametrizing the set 

{(b (Qo-H- J) 6- x )< :teB„) . 

Note that, for i > j we have 

(6 £y r 1 )^ = ((& ei)(ej b' 1 ))^ = {uv T )< , 

where 

U = (n* - n j _ 1 )(6 ei), v T = ( e J &- 1 )(n, - n^i) . 

Thus, 

k 

(b (Qo-H- J) 6- x )<o = ^ (u««a )<o (13-6) 

a=l 

with 

u Q = (Tl ia - n ia _ 1 _ £a _ 1 _ 1 )(& e ia ), v T a = (ef o _ 1 _ Ea _ 1 fT 1 )^ - O^-e^-i) . 
Entries of vectors u a ,v a cannot be arbitrary. First, 

vlu a = ej a b-\U ta - U Ja ^))be ta = ej a e la = . (13-7) 

Next, if e a = 0, i.e. j a = i a -i, then 

Kk = = ■ ( 13 - 8 ) 

Finally, if e Q = 1, i.e. j a — i a -i — 1, then 

wJu Q _l = («a)i a (Wa-l)» a _i-l + («a)i a +l(«a-l)» a _i = (^kjeAWa+l + (^^Jc^ + l^c + ljc + l = • 

(13-9) 

We claim that H13-7fl . (|13-8fl . (|13-9fl are the only restrictions on u a ,v a . We will verify this claim for k = 2. 
The general case follows by an easy induction. 

If ei = 0, we set ui = colfttn, U12, U13, 0, . . . , 0] and vj = [v%\, vf 2 , V13, 0, . . . , 0], where Uu, W13 7^ 
0, W11 7^ 0, V13 S C and U12, U12 G C 4l ~ 2 . Similarly, u 2 = col[0, . . . , 0, U21, U22, U23] and = [0, . . . , 0, V21 = 
u i3' v 22: ^23], where 1121,^23 7^ 0,^23 € C and 1/22,^22 € C n_ll_1 . We assume that conditions i|13-7|) are 
satisfied: vfui — v 2 u 2 = and define 
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The specified entries are consistent with the relation bb~ x — 1 and entries marked by *s are uniquely 
determined by this relation. 

Similarly, if e\ = 1, we set u\ = col[un, u%2, U13, U14, 0, . . . , 0] and vf = [vn, v\ 2 i ^13, W14, 0, . . . , 0], 
where un,«i 3 , u u 7^ 0, v-y-y 7^ 0, ui 3 , v u e C and W12, W12 € C' 1-3 ; and w 2 = col[0, . . . , 0, u 2 \, u 2 2, U23, "24] 
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and uj = [0, . . . ,0, v 2 i,v 22 = -v 2 i^,vj 3 ,v 2i ], where u 2 i,u 22 ,U 2 4, ^ 0,u 2 i ^ 0,d 2 4 G C and u 2 3,«23 G 
C" _ll_2 . Assuming again that v\u\ — v^u 2 = 0, define 



6 = 





Cl 1 


— 1 T 
«12 




Mil 








\ 




( 


Vll 


T 


l>13 


1)14 


* 









1 





U12 

















1 





— Ml2M]~4 


* 












"a 


Ul3 


— 1 T 
-«21 ^23 


«ai 














'■21 


V22 


T 
«23 


V21 













«14 





■U22 















U li 





— U22U 
















1 


M23 




















1 


~U23U 


I 

















U24 


/ 





















U 2i 



and observe that (|13-7|) . I|13-9[l are consistent with bb~ 1 = 1 and entries marked by *s can be uniquely 
determined . 

To conclude the proof, observe that the right hand side of l(13-6[) is invariant under a transformation 
u a — > t a u a , v a — > t~ v a , where t a are arbitrary non-zero parameters. Therefore, we can assume that 

i v a)j a =1 if a = I or e Q _i = 1 (13-10) 

Recall, that, if e a -i — 0, then (v a )j a is given by l|13-8|) . while (u a )j a is determined by the condition l|13-7|l 
for all a. Furthermore, if e a —i = 1, then j a = i Q _i + l and (j a ,ja), (ja + l> ja) an( A (ja + 1, ja + l)-entries 
of the right hand side of i|13-6|) are given by 

(jaja) ■ (U a -l)ja,(v a -l)j a +( u a)j a +l (£l~ 1 ) ) /j' +1 ~ E«=j a +2 ("«)«("<»)» 

(j'a + l)ia) : K*-l)j„+l(«a-l)j a + ( u a)j a +l l 13 " 11 ) 

(ja + 1, ja + 1) : K*-l)j Q +lK*-l)j Q + l - ( M a)j c ,+l ( ^°; ) 1 ^° 1 , 

where we have used lfl3~7j) . ifHHty and (13-9). Note that the entries in (13-10) are the only entries 
in 113-6(1 that depend on (v a —i)j a ,(v a —i)j a and (u a )j a +i- Moreover, (13-10) does not change under a 
transformation 

(U a )j a +1 — > ( u a)j„ + l - t(u a -i)j a + i , (v a -l)j a — > («a-l)j a + t , («a-l)i Q + l — ► («a-l)j a ~ ~ ■ 

This means that we can set 

(u a )j a +i = (u a )i a _ 1 =0 or e Q _! = 1 (13-12) 

Under the normalizations ((13-10(1 . (|13-12|) and restrictions l(13-7|l . l(13-8|l . l(13-9|l . the rest of the parameters 
in HSU, 

(^a)s? \U(x)s i $ ^Q — 1 4" 1 , ■ ■ ■ , ^ Q , tl — 1, . . . , k , 

can be chosen arbitrarily and, on the other hand, these parameters are uniquely determined by the 
right hand side of i|13-6|) . Thus, for Qq satisfying conditions of Theorem 113. II we have found an explicit 
parametrization of Oq by 2n — 2 independent parameters, which completes the proof. Q.E.D. 
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